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ABSTRACT 

Current cosmological observations show a strong signature of the existence of a dark en- 
ergy component with negative pressure. The most obvious candidate for this dark energy is 
the cosmological constant (with the equation of state wx ~ p/p ~ —1), which, however, 
raises several theoretical difficulties. This has led to models for dark energy component which 
evolves with time. We discuss certain questions related to the determination of the nature of 
dark energy component from observations of high redshift supernova. The main results of 
our analysis are: (i) Even if the precise value of wx is known from observations, it is not 
possible to determine the nature of the unknown dark energy source using only kinematical 
and geometrical measurements. We have given explicit examples to show that different types 
of sources can give rise to a given wx ■ (ii) Although the full data set of supernova observa- 
tions (which are currently available) strongly rule out models without dark energy, the high 
(z > 0.25) and low (z < 0.25) redshift data sets, individually, admit decelerating models 
with zero dark energy. Any possible evolution in the absolute magnitude of the supernovae, if 
detected, might allow the decelerating models to be consistent with the data, (iii) We have in- 
troduced two parameters, which can be obtained entirely from theory, to study the sensitivity 
of the luminosity distance on wx- Using these two parameters, we have argued that although 
one can determine the present value of wx accurately from the data, one cannot constrain the 
evolution of wx ■ 



1 INTRODUCTION 

■ One of the equations governing the dynamics of a Friedman uni- 
' verse, (a/a) — — (47rG/3)(p + 3p), implies that the universe 
will accelerate (a > 0) if (p + 3p) < 0. The analysis of high- 
redshift supernova data (Riess et al. 1998; Perlmutter et al. 1999; 
Riess 2000) seems to suggest that a is indeed positive thereby re- 
quiring (p + 3p) < 0. This condition requires at least p or p to 
be negative. Cosmologists have not yet become desperate enough 
to suggest p < 0; but there was remarkably low resistance in 
the community to accepting the existence of a constituent with 
w = (p/p) < —(1/3). This was, to a great extent, facilitated 
by the fact that cosmologists have long since toyed with the ex- 
istence of a non-zero cosmological constant with the stress-tensor 
Tk = PaSI corresponding to wa = —1. In fact, much before the 
arrival of supernova data on the scene, there were indications that 
the universe actually may be populated by a energy density which is 
very smoothly distributed over large scales. A very clear argument 
to this effect was given based on the analysis of APM data in 1990 
(Efstathiou, Sutherland, & Maddox 1990). Later, the first analysis 
of the COBE data in 1992 indicated that the shape of the standard 
cold dark matter (SCDM) spectrum needs to be modified and the 
existence of a smoother energy distribution, like the cosmological 
constant will be required (Padmanabhan & Narasimha 1992; Efs- 
tathiou, Bond, & White 1992). An analysis of a host of observations 
available in 1996 was used to suggest that the data supports the 
existence of a non-zero cosmological constant (Bagla, Padmanab- 
han, & Narlikar 1996). Given this background, it was not surprising 



that the cosmologists were not too shocked by the supernova data 
which became available from 1998 onwards. Currently, there are 
also other observations, like those of the age of the universe, grav- 
itational lensing surveys etc., which strongly suggest the presence 
of a positive cosmological constant [see Padmanabhan (2002a) for 
a comprehensive review and references]. 

There are, however, well known deep theoretical problems 
with the existence of a non-zero cosmological constant A with a 
magnitude of about K(Gh/c?) ~ 10"^^"'. [These are well docu- 
mented in the literature and will not be discussed here; for a re- 
cent review see Padmanabhan (2002a).] This has prompted a host 
of activity in which one looks for a dark energy component in the 
universe with w < —(1/3) which is different from cosmological 
constant [see Padmanabhan (2002a) and references cited therein]. 
If Wx (a) denotes the time dependent equation of state parameter 
for an unknown dark energy component in the universe, then the 
equation of motion d(pxa'^) ~ wxpxda'^ integrates to give 

px{a) = px{ao)expi^-3 J ^[l + wx(a')]| (1) 

If Wx is not identically equal to —1, then the dark energy den- 
sity will evolve with time (even if wx is a constant). This suggests 
the possibility that the cosmological "constant" can be time depen- 
dent thereby alleviating some of the difficulties which arises in the 
case of Wx = — 1- Several such models (usually based on scalar 
fields) have been constructed in the literature, all of which, gener- 
ically have a time dependent wx{a,). The introduction of such an 
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ad-hoc, time dependent function into cosmology, of course, takes 
away much of the credibility; nevertheless, this issue needs to be 
settled ultimately by observations and not by aesthetics. (If the ob- 
servations forces us into a solution which is unaesthetic by current 
standards, we will either conveniently change the standards or live 
with it!) 

In this paper, we discuss certain questions related to the deter- 
mination of the nature of dark energy component from observations 
related to supernova. This issue has been studied in great detail by 
several groups and our contribution in this paper will border on ped- 
agogy and will present a particular point of view. In the next section 
we briefly recall and stress some inherent theoretical degeneracies 
in these analysis. In Section 3 we reanalyze the currently available 
supernova data in order to focus attention on some key elements. 
(These results exist in alternate forms in published literature but 
we beheve our analysis brings out some features clearly.) Based on 
the lessons learnt in this section, we carry out a similar analysis 
for possible future supernova observations and point out what can 
and cannot be achieved. We have intentionally kept the data analy- 
sis at a fairly simple level (indicated in the title by 'a theoretician's 
analysis'). We subscribe to the point of view that any result which 
cannot be revealed by a simple analysis of data, but arises through 
a more complex statistical procedure, is inherently suspect and a 
conclusion as important as the existence of dark energy with p < 
should pass such a test. 



2 THEORETICAL DEGENERACIES IN THE 
FRIEDMANN MODEL 

We begin by recalling and stressing some inherent limitations 
which exists in all attempts to probe the universe through geometri- 
cal measures. The assumption of isotropy and homogeneity implies 
that the large scale geometry of the universe can be described by a 
metric of the form 



ds^ = — o^(t)dx^ 



(2) 



where dx'^ denotes the line element of the three-dimensional space 
in comoving coordinates. In any range of time during which a(t) is 
a monotoiric function of t, one can use a (or equivalently, redshift 
z = {a/ao)~^ — 1) as a time coordinate. The metric is then 
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where H{z) = d/a is the Hubble parameter. 

This equation allows us to draw the following conclusion: The 
only non-trivial metric function in a Friedmann universe is the func- 
tion H{z) (besides the curvature of the spatial part of the metric). 
Hence, any kind of observation based on geometry of spacetime, 
however complex it may be, will not allow us to determine anything 
other than this single function H{z). Since Friedmann equations 
relate H^{z) to the total energy density in the universe (assuming 
that the curvature of the spatial part of the metric is known from in- 
dependent observations or fixed by some theoretical prejudice), the 
best we can do from any geometrical observation is to determine 
the total energy density of the universe at any given «. It is not pos- 
sible to determine the energy densities of individual components 
from any geometrical observation. 



More explicitly, when several non-interacting sources are 
present in the universe, the Friedmann equations give 
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where Q,a is the density parameter for the a-th component (Uke ra- 
diation, matter, cosmological constant etc.) and Wa{z) is the corre- 
sponding equation of state parameter. For example, non-relativistic 
matter {Q.m) has Wm = 0, while a non-evolving cosmological con- 
stant (Qa) has wa = — 1. From the above equation, it is clear 
that one needs to know the Sla and Wc, (z) of the N — 1 compo- 
nents [along with geometrical measurements giving H{z)] in order 
to determine the contribution of the iV-th component. If we take 
Qk = corresponding to the flat universe, and assume that there is 
only dust-like matter and dark energy present in the universe, then 
we need to know Urn of non-relativistic dust-like matter and H{z) 
to get a handle on fix of the dark energy. Unfortunately, even after 
making these assumptions, we are plagued by the uncertainties in 
Qrn which is currently estimated to be anywhere between 0.2 and 
0.35. As has been noted several times in the literature (and as we 
shall emphasize in Section 4), this is a fairly strong degeneracy. 

There is another — and from a theoretical point of view more 
serious — degeneracy which seems to have been inadequately 
stressed in the literature. Let us assume that the universe is made 
of two components: Pkn(a), which is known from independent ob- 
servations and a component px (a) which is not known. From the 
Friedmann equation, it follows that 
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Taking a derivative of Inpx(a) and using (1), it is easy to obtain 
the relation 
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3 din a 

If geometrical observations of the universe give us Hia) and other 
observations give us /9kn(ci) then one can determine nkn(a) and 
thus TOx(ffl). This is the best uniform background cosmology can 
do for us. As far as cosmological model building is concerned, one 
is usually interested in the features of the dark energy "fluid", and 
it is usually sufficient to know wx (a) and the speed of sound of 
the fluid. However, theoretical physicists would like to know some- 
thing more about the nature of the dark energy "fluid" - in partic- 
ular, they would find it useful to determine the Lagrangian for the 
fluid. (For example, in the case of inflation, a considerable amotmt 
of effort is being spent to understand the form of the potential for 
the scalar field.) In this paper, we shall take the viewpoint of a the- 
oretical physicist, where the nature of the dark energy is defined by 
its Lagrangian. 

Now, if a geometrical observation in future suggest that 
wx{a) is not identically equal to —1, then the question arises as 
to what is the nature of this dark energy field. It seems virtually 
impossible that its nature can be determined from laboratory exper- 
iments (unlike, for example, WIMPS which might constitute dark 
matter). However, knowing wx (a) is grossly inadequate for deter- 
mining the physical nature of dark energy. For example, even if one 
makes another gigantic leap of faith and assumes that the dark en- 
ergy arises from a scalar field, it is possible to come up with scalar 
field Lagrangians of different forms leading to same wx (a) [see 



A theoretician's analysis of the supernova data 3 



Quintessence Tachyon 




Ho ^ 

Figure 1. The scalar field potentials for the quintessence (left frame) and tachyonic (right frame) fields which are required to produce the dark energy equation 
of state wx{0') = "Jo — wi(a — 1). The potentials are plotted for the parameter values Q-m = 0.3, ilx = 1 ~ = 0.7, wq = 0.5, wi = —0.1 using 
equations (10) - (13). Since the equations give the potential V{(f>) in the parametric form [^(z), (f>{z)], each point on the V — </> curve can be labelled by the 
corresponding value of redshift, z. We have pointed some particular redshifts in both the curves. It is clear that the scalar field starts from a large value of the 
potential, roUs down as time progresses and the density contributed by the potential reaches a value ~ Clx aXz = 0. 



Padmanabhan (2002b)]. To illustrate this point, we shall discuss 
two possibilities : 

1 , 



Lquin = -di4>d'<t> - V{4>), Ltach = -V{4>)^/l^^a4d^ (7) 

Both these Lagrangians involve one arbitrary function ¥{4)). The 
first one, Lquin, which is a natural generalization of the Lagrangian 
for a non-relativistic particle [L — /2 — V{q)}, is usually called 
quintessence. Similarly, the second one, Ltach, is a generalization 
of the Lagrangian for a relativistic particle [L = — my'l — q^], 
and is usually called a tachyonic Lagrangian. [It arises naturally in 
the string theoretical context; see Sen (2002a); Sen (2002b); Sen 
(2002c). For a discussion of the cosmological aspect of this La- 
grangian, see Padmanabhan & Choudhury (2002) and references 
therein.] When these Lagrangians act as sources in Friedmann uni- 
verse, they are characterized by density and equation of state pa- 
rameters given by 
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respectively. 

Since both the Lagrangians have one undetermined function 
V{4>), it is possible to choose this function in order to produce a 
given wx{a,) [or equivalently, H{a)]. For a flat universe, one can 
determine the function V{(p) for the quintessence model from the 
implicit relations (ElUs & Madsen 1991) 
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Similarly, in the case of tachyonic scalar field, the potential function 
is given by (Padmanabhan 2002b) 
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,3[1 -nu„(a)] 3 H{a) 

We can present an explicit example of a situation where wx 
is assumed to be known of the simple form 

z 



wx{z) = Wo — wi(a — 1) = Wo + wi 



1 + z' 



(14) 



where wo measures the current value of the parameter and — wi 
gives its rate of change (with respect to the scale factor) at the 
present epoch. In addition to simplicity, this parametrization has 
the advantage of giving finite wx in the entire range < « < oo. 
We shall use this simple parametrization later to discuss the possi- 
bility of constraining wx from observations. The Hubble parameter 
in this case is given by 



+ {1-Q.m){l + Z) 
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(15) 
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This equation of state parameter can be obtained (for small val- 
ues of I wi I) from at least two scalar field Lagrangians - those for 
the quintessence and tachyoiuc fields. The corresponding poten- 
tials V{(t>) can be determined from equations (11) and (12). Un- 
fortunately, the relations cannot be written in a closed form (un- 
less w\ = 0) in this case. The form of the potential for both the 
cases can be calculated numerically for given values of wo and 
wi, which are plotted in Figure 1. The left frame shows V{4i) for 
the quintessence field while the right frame shows that for the and 
tachyonic field. The potentials arc plotted for the parameter val- 
ues Q,jn = 0.3, fix = 1 — = 0.7, Wo = 0.5, Wl = —0.1 
using equations (10) - (13). Since the equations give the potential 
V{4>) in the parametric form [1^(2;), </'(z)], each point on the V — 4> 
curve can be labelled by the corresponding value of redshift, z. We 
have pointed some particular redshifts in both the curves. It is clear 
that the scalar field starts from a large value of the potential, rolls 
down as time progresses and the density contributed by the poten- 
tial reaches a value ~ Q.x at « = 0. 

This discussion shows that even when wx{a) or H{a) is 
known, it is not possible to proceed further and determine the nature 
of the unknown dark energy source using only kinematical and geo- 
metrical measurements. This constitutes a serious theoretical prob- 
lem in cosmology if observations suggest that wx is not identically 
equal to —1 (Huterer & Turner 1999). 



3 WARM UP: CURRENT SUPERNOVA DATA AND 
THEIR ANALYSIS 

We shall next reanalyze the currently available supernova data in 
order to stress certain features which are inherent in such an analy- 
sis. This may be considered a warm up exercise for the next section. 

The observations directly measure the apparent magnitude rn 
of a supernova and its redshift z. The apparent magnitude m is 
related to the luminosity distance of the supernova through 



m{z) = M + 5 log 



dL{z) 



1 Mpc 



+ 25, 



(16) 



where M is the absolute magnitude (which is believed to be con- 
stant for all supemovae of Type-la - this is what is called the 
"standard candle hypothesis"). It is convenient to work with a di- 
mensionless quantity (called the "Hubble-constant-free" luminos- 
ity distance) 



Q{z) = ^dL{z) 
which gives 

m(0) = M -|-51ogio(5(«), 
where 
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M-51ogi„/i + 42.38.(19) 



Any model of cosmology will predict Q{z) with some unde- 
termined parameters (say, for example, Q.m,^A)- These parame- 
ters, along with M (which is itself a combination of the absolute 
magnitude of the Type-la supernova and the Hubble constant), are 
obtained by comparing with observations. According to the stan- 
dard hot big bang model of cosmology [see e.g., Padmanabhan 
(2002c)], 
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Figure 2. Comparison between various flat models and the observational 
data. The observational data points, shown with error-bars, are obtained 
from Perhnutter et al. (1999) and Riess et al. (2001). In order to determine 
Q(z) for these data points, we have used equation (18) with the best-fit 
value M = 23.95. 
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The data in this paper is based on the redshift-magnitude re- 
lation of 54 Type-la supemovae [excluding 6 outHers from the full 
sample of 60, which consists of 18 low-redshift points from the 
Calan/Tololo supernova survey (Hamuy et al. 1996) and 42 points 
from the Supernova Cosmology Project (Perlmutter et al. 1999)] 
and that of supernova 1997ff at « = 1.755 (Riess et al. 2001). The 
values of m used have been already corrected for the supernova 
light curve width-luminosity relation, galactic extinction and pos- 
sible K-correction. In addition, the magnitude for supernova 1997ff 
has been corrected for lensing effects (Benitez et al. 2002). Given a 
cosmological model with some free parameters, we have obtained 
the best-fit parameter values and the corresponding covariance ma- 
trix using the Levenberg-Marquardt method (Press et al. 1992). 
The details of this fitting procedure, although quite well-known, 
are given in Appendix A for completeness. 

Let us start with flat models with Q.m -h Ha = 1; fife = 0. 
which are currently favoured strongly by CMBR data (for recent 
results, see Sievers et al. 2002). The simple analysis mentioned 
above gives a best-fit value of Q,m (after marginalizing over M) 
to be 0.31 ± 0.08. The best-fit M. (after marginalizing over Q,m) 
is found to be 23.95 ± 0.05 (all the errors quoted in this paper are 
la). The comparison between three flat models and the observa- 
tional data is shown in in Figure 2. In order to determine Q{z) for 
the observational data points, we have used equation (18) with the 
best-fit value M = 23.95. 

Although the best-fit analysis shows that the data favour 
strongly for a positive non-zero f^A which, in turn, implies the pres- 
ence of an accelerating universe, the same conclusion is not visually 
obvious from Figure 2. In order to see that the data favours mod- 
els with non-zero Qa, one usually plots the supernova magnitude 
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Figure 3. The observed supernova data points in tlie a — a plane 
for flat models. The procedure for obtaining the data points and the 
corresponding error-bars are described in the text. The sohd curves, 
from bottom to top, are for flat cosmological models with Qrn = 
0.00, 0.16, 0.32, 0.48, 0.64, 0.80, 1.00 respectively. 



with respect to a fiducial best-fit model [see, for example, Perlmut- 
ter et al. (1999)] - however, to see the presence of an accelerating 
phase, it is more convenient to display the data as the phase portrait 
of the universe in the a — a plane. (The procedure for doing this 
is described in Appendix B; as far as we are aware, this has not 
been done in hterature before.) The data points, with error-bars, in 
the a — a plane are shown in Figure 3. The solid curves plotted in 
Figure 3 correspond to theoretical flat models with different Q.m- 
The data points show a clear sign of an accelerating universe at low 
redshifts. Hence, in principle, one should be able to rule out non- 
accelerating models using only the low redshift data. However, as 
is clear from Figure 3, it is not possible to rule out any of the cos- 
mological models using low redshift {z < 0.25) data because of 
large error-bars. (The fact that low redshift data cannot rule out any 
of the cosmological models can be seen in Figure 2 also.) On the 
other hand, high redshift data alone cannot be used to establish the 
existence of a cosmological constant. For example, one can use the 
freedom in the value of M to shift the data points vertically, and 
make them consistent with a decelerating model (Clm = 1, topmost 
curve). It is the interplay between both the high and low redshift 
supernova data which leads to a clear indication of an accelerating 
phase. 

All the above conclusions can be made more quantitative by 
studying at the confidence ellipses in the Qm — M plane, shown 
in Figure 4. The method of drawing these ellipses is outlined in 
Appendix A. In the left panel, we have plotted the confidence re- 
gions using the full data set of 55 points. It is obvious that most of 
the probability is concentrated around the best-fit value. The confi- 
dence contours in the middle and right panel are obtained by repeat- 
ing the best-fit analysis for the low redshift data set (z < 0.25) and 
high redshift data set (z > 0.25), respectively. It is clear from the 
middle panel that, although the low redshift data can constrain M 
very well, it is unable to constrain Qrn- On the other hand, the high 



redshift data (right panel) is able to constrain neither M nor Qm ■ 
In particular, the decelerating model (Qm = 1) is quite consistent 
with both the low and high redshift data sets when they are treated 
separately. One needs to combine the low and high redshift data to 
constrain Qm ~ because of the angular slant of the ellipses, a best- 
fit region around Qm = 0.31 is isolated, as can be seen from the 
left panel. [This analysis indirectly stresses the importance of any 
evolutionary effects. If, for example, supernova at 2 > 0.25 and su- 
pernova at z < 0.25 have different absolute luminosities because 
of some unknown effect, then the entire data set can be made con- 
sistent with Qm = 1, = model. It certainly appears ad-hoc; 
but one should compare the ad-hocness in any of these assump- 
tions with the ad-hocness in introducing a dimensionless constant 
A{Gh/c^) ~ 10^^^^ in the physical system to explain the cosmo- 
logical observations (Padmanabhan 2002a)!] 

We now generalize some of the results for non-flat cosmolo- 
gies. In this case, we have three free parameters, namely, Qm, Qa 
and A4. The confidence region ellipses in the Qm-QA plane can 
be drawn in two ways - (i) by keeping Al as a free parameter and 
then marginalize the joint probability distribution over A4 or, (ii) 
by fixing M to some constant value (and deal with only two free 
parameters, Qm and Qa). 

The confidence region ellipses in the Qm-QA plane are shown 
in Figure 5 for the full data set. In panel (a) the confidence re- 
gions are obtained by marginalizing over M. The best-fit values 
are found to be Qm = 0.67 ± 0.25, Qa = 1.24 ± 0.34 and 
A4 = 23.92±0.05, as indicated in the panel. For panels (b - d), the 
confidence contours are obtained by fixing to a constant value 
rather than marginalizing over this parameter. The three frames cor- 
respond to the best-fit mean value and two values in the wings of 
la from the mean, respectively. The fixed value of A4 is indicated 
in the panel. 

The main conclusions we can draw from this figure are: (i) The 
results do not change significantly whether we marginalize over M. 
or whether we use the best-fit value. One can see that the best-fit 
values of Qm and Qa do not change at all, it is only the the spread 
which decreases slightly when we fix Af to a particular value rather 
than marginalize it. The two methods give such sinular results be- 
cause the probabiUty of M is sharply peaked (the spread is ~ 1 
per cent of the mean value), (ii) The results exclude the SCDM 
model (Qm = 1 , Qa ~ 0) at a high level of significance, taking 
into account the maximum uncertainty in A4. The importance of 
this exercise will be evident later when we consider high and low 
redshift data points separately, (iii) The slanted shape of the proba- 
bility ellipses show that a particular linear combination of Qm and 
Qa is selected out by these observations (which, in this case, turns 
out to be 0.82f2,n — 0.57f2A). This feature, of course, has nothing 
to do with supernova data and arises purely because the luminos- 
ity distance Q depends strongly on a particular linear combination 
of Qm and Qa - This point is illustrated by plotting the contour of 
constant luminosity distance, Q{z — 0.63) = constant. The coin- 
cidence of this line (which roughly corresponds to Q at a redshift in 
the middle of the data) with the probability ellipses indicates that 
it is the dependence of the luminosity distance on cosmological 
parameters which essentially determines the nature of this result. 
This fact is known in the hterature [see e.g., Goobar & Perlmutter 
(1995)], but we have not seen the actual data presented in this form. 
We shall use this result in the next section to discuss the possibility 
of determining the equation of state parameter for the evolving dark 
energy component. One should note that a complicated likelihood 
analysis may give confidence contours of a different shape (they 
are ellipses only because of the assumption of normal distribution 
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Figure 4. Confidence region elHp.ses in the f2„i — M plane for flat models witli non-relativi.stic matter and a cosmological constant. Tlie ellipses corresponding 
to the 68, 90 and 99 per cent confidence regions are shown. In the left panel, all the 55 data points from Perlmutter et al. (1999) and Reiss et al. (2001) are 
used. In the middle panel, data points with z < 0.25 are used, while in the light panel, we have used data points with z > 0.25. We have indicated the best-fit 
values of and JV[ (with \a en'ors). 
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Figure 5. Confidence region ellipses in the ilm — f^A plane for models with non-relativistic matter and a cosmological constant. The ellipses corresponding 
to the 68, 90 and 99 per cent confidence regions are shown, (a) The confidence regions are obtained after marginalizing over A4. The best-fit values of f2„i 
and M are indicated (with Icr eiTors). (b - d) The confidence regions are obtained after fixing J\A to its best-fit mean value and the to the values in the wings 
of Icr limit respectively. The fixed value of A4 is indicated in the panel. The dashed line corresponds to the flat model (Hm + = 1). We have indicated the 
value of for the best-fit parameters in all the panels. The unbroken slanted line corresponds to the contour of constant luminosity distance, Q{z = 0.63) = 
constant. 
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of errors here) - however, the shapes of those contours are bound 
show the degeneracy in Q.m and (arising purely from theory) 
mentioned above. We should also mention that the shapes of the 
contours of constant Q at high redshifts are quite different from 
those at low redshifts - hence more high redshift data can be useful 
in breaking the degeneracy between Q,m and Oa. 

As discussed in the third point above, the luminosity distance 
is sensitive only to a particular linear combination ofQ,m and JIa, 
which is illustrated in Figure 6. In this figure, Q,m and Qa are 
treated as free parameters (within 2(j bounds from the best-fit val- 
ues mentioned above), but the combination Q.&2Q,m — 0.57fiA is 
held fixed. It turns out that Q{z) is not very sensitive to individual 
values of Q,m and JIa at low redshifts when 0.8251^ — 0.57f2A is 
in the range —0.14 ± 0.11. This is clear from Figure 6 in which 
a wide variety of cosmological models are plotted along with the 
data for a constant value of the above combination. Though some 
of the models have unacceptable values of flm and Q.h (ruled out 
by observations of CMBR and structure formation), the supernova 
measurements alone cannot rule them out. Essentially, the data at 
z < 1 is degenerate on the linear combination 0.820^ — 0.57f2A. 
Our analysis of the supernova data shows that the best-fit value for 
this combination is 0.82fim - 0.57Qa = -0.14 ± 0.11. 

Finally, we comment on the interplay between high and low 
redshift data for non-flat models. Just as in the case of the flat 
models, we divide the full data set into low {z < 0.25) and high 
(z > 0.25) redshift subsets, and repeated the best-fit analysis. The 
resulting confidence contours are shown in Figure 7. In panel (a), 
the confidence contours are plotted for the high redshift data set af- 
ter marginalizing over A4. As we stressed before, the SCDM model 
cannot be ruled out using high redshift data alone. In panels (b - d), 
we show the corresponding results in which the values of A4 are 
fixed (rather than marginalizing over M). Comparing the frames 
(a - d) of Figure 7 with the frames (a - d) of Figure 5 where the 
full data set was used, we can draw the following conclusions: (i) 
The best-fit value for M is now 24.05 ± 0.38; the la error has 
now gone up by nearly an order of magnitude compared to the case 
where the full data was used. Because of this spread, the confidence 
contours are sensitive to the value of A4 one uses, unlike the situ- 
ation where all the data points were used, (ii) Our conclusions will 
now depend on the exact value of M used. For the best-fit mean 
value and the lower end of A4 , the high redshift data can rule out 
the SCDM model [see frames (b) and (d)]. But for the higher end of 
the allowed la range of TV), we cannot exclude the SCDM model 
[see frame (c)]. This essentially shows that it is difficult to rule out 
models when there are large tmcertainties in M. Finally, we plot 
the contours for the low redshift data set after marginalizing over 
A4 in panel (e). It is clear that the low redshift data cannot be used 
to discriminate between cosmological models effectively. This is 
because Q{z) at low redshifts is only very weakly dependent on 
the cosmological parameters. So, even though the acceleration of 
the universe is a low redshift phenomenon, we cannot reliably de- 
termine it using low redshift data alone. 



4 CONSTRAINTS ON EVOLVING DARK ENERGY 

As we have seen in the previous two sections, supernova obser- 
vations suggest the existence of a non-zero Oa. However, this very 
existence of JIa raises serious theoretical problems [see Padmanab- 
han (2002a) for a detailed discussion]. These difficulties have led 
people to consider the possibility that the dark energy is not just 
a constant - but is evolving with time. The evolution of the dark 
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Figure 6. The luminosity distance for a class of models with Sim varying 
between 0.67±0.50 and SIa between 1.24±0.74, but with a constant value 
of the combination 0.820^ — 0.57f2A- It is clear that when the combina- 
tion 0.82Qm — 0.57f2A is fixed, low redshift observations cannot distin- 
guish between the different models even if and Qa vary significantly. 



energy component can be parametrized by its equation of state pa- 
rameter wx(z), the evolution of the density px being given by 
equation (1). In this section, we shall examine the possibility of 
constraining wx («) by comparing theoretical models with super- 
nova observations. 

One simple, phenomenological, procedure for comparing ob- 
servations with theory is to parametrize the function wx{z) in 
some suitable form and determine a finite set of parameters in this 
function using the observations. Theoretical models can then be re- 
duced to a finite set of parameters which can be determined from 
observations. To illustrate this approach, let us assume that, in the 
flat universe, wx {z) is completely determined by its current value 
(wo) and its rate of change with respect to the scale factor at the 
present epoch (— wi) and is given by the simple form as in equa- 
tion (14). 

Since the cosmological constant corresponds to the case wq = 

— l,w\ = 0, it is interesting to see the constraints on wa from 
supernova data even if we assume w\ = 0. In such a study, it should 
be noted that (i) the acceleration of the universe requires wq < 

— 1/3, and (ii) all values of of wq other that wo = —1 leads to 
a dark energy density which evolves as (1 + z)"^'^^""-*. A simple 
best-fit analysis shows that for a flat model with f2rn = 0.31 and 
M. = 23.95 (the best-fit parameters for flat models, obtained in the 
previous section), the best-fit value of wo is —1.01 ± 0.12 (which 
is nothing but the conventional cosmological constant). The data 
clearly rules out models with wo > —1/3 at a high significance 
level, thereby supporting the existence of a dark energy component 
with negative pressiu'e. 

One can extend the analysis to find the constraints in the 
Wo — wi plane. [There is extensive literature on determining dark 
energy parameters from the current and proposed future super- 
nova observations, e.g., see Maor et al. (2002); Weller & Albrecht 
(2002); Gerke & Efstathiou (2002) and references therein.] As be- 
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Figure 7. Confidence region ellipses in the Q,„ — Ha plane for models with non-relativistic matter and a cosmological constant. The ellipses corresponding 
to the 68, 90 and 99 per cent confidence regions are shown, (a) The confidence regions are obtained for high redshift points (z > 0.25), after marginalizing 
over M. The best-fit value (with Itr error) of JVl is indicated, (b - d) The confidence regions are obtained for high redshift points (z > 0.25), after fixing A4 
to its best-fit value, the upper limit and the lower limit respectively. The fixed value of A4 is indicated in the panel, (e) The confidence regions are obtained for 
low redshift points {z < 0.25), after marginalizing over M. The best-fit value (with \a error) of M is indicated. 




Figure 8. Confidence region ellipses in the ioq — fi plane for flat models with Q,m = 0.21, 0.31 and 0.41 respectively, as indicated in the frames. The value 
of M is chosen to be the best-fit value, which is also indicated. The ellipses corresponding to the 68, 90 and 99 per cent confidence regions are shown. The 
square point denotes the equation of state for a universe with a non-evolving dark energy component (the cosmological constant). The unbroken slanted line 
corresponds to the contour of constant luminosity distance, Q(z = 0.6) = constant. 
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'Wo 

Figure 9. Sensitivity of Q{z) to tlie parameters wq and wi . Tlie curves correspond to tlie constant values for tlie fractional change in Q for unit change in wq 
(top frames) and uii (bottom frames) for two redshifts 2 = 1 (left frames) and 2 = 3 (right frames). The value of Q for a particular curve is indicated. A flat 
cosmological model with Sim = 0.31 has been used. The shaded bands across the frames correspond to the regions in which —1 < wx{z) < 0. 



fore, we assume a flat universe with a fixed value of Sim in the 
range (0.21 - 0.41) and M is fixed to the corresponding best-fit 
value. The confidence contours for the three models are shown 
in Figure 8. The square point denotes the equation of state for a 
universe with a non-evolving dark energy component (the cosmo- 
logical constant). The main points revealed by this figure are: (i) 
The equation of state corresponding to the cosmological constant 
is within the la contour for 0.21 < Qrn < 0.41. (ii) Despite the 
uncertainties in wi, models with wq > —1/3 are ruled out at high 
significance level for Q,m < 0.4. (One can accommodate models 
with Wo > —1/3 only for Q,rn > 0.4 for very high negative values 
of wi.) (iii) The shape of the confidence contours clearly indicates 
that the data is not as sensitive to wi as compared to wq. Since the 
supernova observations essentially measure Q{z), this is a clear in- 
dication that Q{z) is comparatively insensitive to wi. To see this, 
note that in Figure 8, the position of the slanted line (which cor- 
responds to Q{z = 0.6) = constant) exactly coincides with the 
contour ellipses. We mention again that the shapes of the confi- 
dence contours may not remain as simple ellipses if one carries out 
a complicated likelihood analysis - however, as our conclusions are 
based on the sensitivity of Q{z) on wq and lui, they are likely to re- 
main the same for any analysis. Since the determination of wq and 
wi is of considerable importance, we shall provide an illustration 
of the relative sensitivities of geometrical features to wq and wi. 

As we have stressed in the previous section and also above, 
the supernova observations essentially measure Q{z), and hence 
the accuracy in the determination of wq and wi from (both present 
and planned future) supernova observations will crucially depend 



on how sensitive Q is to the changes in wq and wi [which is clear 
from the coincidence of the Q — constant line with the probability 
ellipses in Figure 8]. A good measure of the sensitivity is provided 
by the two parameters 

A(z;wo,wi) = -^ln[Q{z;wo,wi)], 
awo 

B{z;wo,wi) = -^^ln[Q{z;wo,wi)]. (22) 
awi 

Like Q{z; wo, Wi), the parameters A and B can be computed for 
a particular cosmological model in a straightforward manner. At 
any given redshift z, we can plot contours of constant A and B in 
the Wo — wi plane. Figure 9 shows the result of such analysis for 
flat models with Q,,n = 0.31. The two frames on the left are for 
z — 1, while those on the right are for jz = 3. The top frames 
show the contours of constant A while the bottom frames show 
the contours of constant B. From the definitions (22) it is clear 
that A and B can be interpreted as the fractional change in Q for 
unit change in wo and wi, respectively. For example, along the line 
marked A = 0.2 (in the top left frame), Q will change by 20 per 
cent for unit change in wq. It is clear from the two top frames that 
for most of the interesting region in the wo — Wi plane, changing 
Wo by unity changes Q by about 10 per cent or more. Comparison 
of the two cases z — 1 and z = 3 (the two top frames) shows 
that the sensitivity is higher at high redshift, as to be expected. The 
shaded bands across the frames correspond to the regions in which 
— 1 < wx{z) < 0, which is of primary interest in constraining 
dark energy with negative pressure. One concludes from the above 
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discussion that determining wq from Q fairly accurately will not be 
too daunting a task. 

The situation, however, is quite different for mi as illustrated 
in the two bottom frames. For the same region of the wo — w\ plane, 
Q changes only by a few per cent when lui changes by unity. This 
means that Q is much less sensitive to w\ than to wq. It is going 
to be significantly more difficult to determine a value for tiii from 
observations of Q (like supernova observations) in the near future. 
Comparison of the two cases z = 1 and z = 3 again shows that the 
sensitivity is somewhat better at high redshifts, but only marginally. 

In the above analysis, we have treated to be a constant. 
The situation is made worse by the fact that Q also depends on Q.m ■ 
If the variation of Q,m mimics that of wq or lui , then one also needs 
to determine the sensitivity of Q to Q,m - Figure 10 shows contours 
of constant Q in the Q,m — wi (top frames) and the ilm — wq (bot- 
tom frames) planes at two redshifts z = 1 (left frames) and z — 3 
(right frames). The two top frames show that if one varies the value 
of Qrn in the allowed range, say, (0.2 - 0.4), one can move along 
the curve of constant Q and induce a fairly large variation (~ 1 - 2) 
in wi. In other words, large changes in wi can be easily compen- 
sated by small changes in Sim while maintaining the same value 
for Q at a given redshift. This shows that the uncertainty in Qm 
introduces further difficulties in determining wi accurately from 
measurements of Q. The two bottom frames show that the situation 
is better for wq. The curves are much less steep and hence vary- 
ing Qrn does not induce large variations in wq. We once again are 
led to the conclusion that unambiguous determination of wi from 
data will be quite difficult. It appears that supernova observations 



may not be of great help in ruling out (non-evolving) cosmological 
constant as the major dark energy component. 



5 DISCUSSION 

The supernova data on the whole, rules out non-accelerating mod- 
els with very high significance level, as noted and stressed by vari- 
ous authors (Riess et al. 1998; Perlmutter et al. 1999; Riess 2000). 
It is thus more or less conclusive that we have to look for some 
form of matter which has negative pressure. However, it is interest- 
ing to note that if we divide the data set into high and low redshift 
subsets, neither of the subsets are able to rule out the decelerating 
models ~ it is only the interplay between the high and low redshift 
data which implies the result quoted above. This analysis indirectly 
stresses the importance of any evolutionary effects. Any possible 
evolution in the absolute magnitude of the supemovae, if detected, 
might allow the decelerating models to be consistent with the data. 
However, our understanding of the complicated physical processes 
in the supernovae is definitely not well enough to discuss the possi- 
bility of detecting supernova evolution in near future. Although the 
intrinsic evolution of the SN is an interesting problem, addressing 
these issues is beyond the scope of this paper. 

The supernova data available today is not enough to fix the 
parameters Qrn and SIa independently. Instead, a particular com- 
bination of the two parameters is chosen out by the data. This de- 
generacy has to do with the dependence of the luminosity distance 
on these parameters. Since supernova observations essentially mea- 
sure the luminosity distance Q{z), the accuracy in the determina- 
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tion of various parameters from supernova observations will cru- 
cially depend on how sensitive Q is to the changes in those param- 
eters. This analysis can be done entirely from theory [using Fisher 
information matrix; see for example Efstathiou (1999)] and, in prin- 
ciple, can provide nice handle on which combinations of parame- 
ters are expected to be constrained from the supernova data. 

The key issue regarding dark energy is whether it is a con- 
stant or whether it is evolving with time. In particular, one is in- 
terested to determine the evolution of the equation of state, wx 
of the dark energy component. To address this question, we have 
taken a simple phenomenological model for wx having the form 
wx{a) = Wo — wi(o — 1). The sensitivity of the luminosity dis- 
tance (and hence, the supernova data) on wo and wi can be mea- 
sured by the two parameters (A and B) introduced in the paper, 
which are essentially the fractional changes in Q{z) for unit change 
in Wo and wi, respectively. The parameters A and B can be ob- 
tained entirely from theory and the conclusions drawn from them 
are valid for current as well as future supernova observations. We 
find that Q{z) is quite sensitive to wq - hence, one can constrain 
the current value of wx quite well. However, Q(z) is compara- 
tively insensitive to wi, thus determining the evolution of wx will 
be a difficult task. The situation is further worsened when we take 
the uncertainties in Qm into account. Our results are consistent with 
other statistical analyses done by Astier (2000), Maor, Brustein, & 
Steinhardt (2001), Weller & Albrecht (2001), Maor et al. (2002), 
WeUer & Albrecht (2002). They have shown that it is possible to 
constrain the present value of wx with future missions like "Su- 
perNovae Acceleration Probe" (SNAP); however, constraining the 
evolution of wx will not be easy using SNAP, particularly when 
Qrn is not known to a high accuracy. 
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APPENDIX A: DATA ANALYSIS 

In this section, we shall outline the method used for obtain- 
ing the best-fit (theoretical) parameters from observational data, 
and draw the corresponding confidence contours. For simplicity, 
we shall assume that the measurement errors are normally dis- 
tributed. Suppose we have M observational data points denoted 
by mobs(ai);i = 1,...,M with corresponding errors <Jm,{zi). 
For a given theoretical model mth(z; Ca) with v free parameters 
Ca; a = one can construct the quantity 

M 



(Al) 



The best-fit parameters Cq 
quantity 



dc 



= 



are obtained by minimizing the above 



(A2) 



using the Levenberg-Marquardt method (Press et al. 1992). One 
can show that the quantity Ax^ = X^(ca) — X^(ca) follows a 
chi-square distribution with i/ degrees of freedom. 

To obtain the confidence intervals, let us first define the curva- 
ture matrix 



Aa3 = 



a2 2 

d X 



dCadC/} 



(A3) 



The covariance matrix Cap is simply the inverse of the curvature 
matrix. The probability distribution of the parameters is given by 



V{ca) = const X exp 



2 X/ - Ca)^a/3(C/3 



(A4) 



Suppose we are interested in the confidence regions for a sub- 
set of, say, v' parameters (y' < v). In that case, the the regions 
are obtained by simply marginalizing the probability distribution 
over the rest v — v' parameters [which is equivalent to integrating 
the probability distribution (A4) over the rest v — v' parameters]. 
This can be done simply by taking the full v y,v covariance matrix 
Ca/3 and copying the intersection of the v' rows and columns cor- 
responding to the parameters of interest into a z/' x i^' matrix C^^. 
The inverse of this will give the corresponding curvature matrix 

(Press et al. 1992). The marginalized probability distribution 
will be simply 



V [ca) = const X exp 



(A5) 



Next, one needs to find the quantity Axp, such that the proba- 
bility of a chi-square variable with v' degrees of freedom being less 
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Figure Bl. The difference between the actual data points and the fit (Bl) to 
the data. 



than is p, where p is the desired confidence limit (e.g., 0.68 or 
0.95). This can be obtained by solving the equation 

P ■■ 



\ r 

2-'/2r(z/'/2) 



(A6) 



where T{z) is the gamma function. For a given p and a set of v' 
parameters, the equation for the boundary of the confidence region 
in the i^'-dimensional space is given by the equation 



For flat models, it straightforward to obtain the Hubble parameter 
from Q(z). In particular, we are interested in the quantity 

-1 



{1 + z 



d 



l + z 



(B5) 



which will enable us to plot the data points in the a — a plane. 
The determination of the corresponding error-bars is a non-trivial 
exercise. In this paper, the errors are obtained from the relation 



+ 



d \Q{z){l + eQ{z)) 



1+2 



,(B6) 



where and eg denote the fractional error in d and Q, respec- 
tively. From equation (18), we have ^q^z) = 0.2 hi 10 a-,n{z), 
where am{z) is the total uncertainty in the observed magnitude. 
Since there is no systematic evolution of the observed (Tm, (z) (see 
Figure Bl), we have, to the lowest approximation, the fractional 
error in d 

taiz) = eqiz) =0.2\nWam{z). (Bl) 

Thus, one can plot the quantity a{z) for the observational data with 

error-bars given by the above equation. 

Note that the analysis assumes that the observed errors in the 
measurement of z (or, equivalently a) are negligible. Strictly speak- 
ing, the errors in d calculated above are strict lower limits - they can 
be slightly higher because of the errors arising from (i) the fitting 
function and (ii) any systematic evolution of the observed am{z). 
However, these effects are unlikely to affect any of the conclusions 
we have drawn in this paper. 



^Xp = (Ca - Cc)A'^g{Cf} - Cf}) 



(A7) 



APPENDIX B: DETERMINATION OF THE HUBBLE 

PARAMETER H{z) FROM SUPERNOVA OBSERVATIONS 

The observational data used in this paper can be fitted by the func- 
tion of simple form 

z{l + 022) 



met (2) = ai + 51ogio 



(Bl) 



1 -I- asz 
where the parameters are 

ai = 23.95 ± 0.05; a2 = 2.00 ± 1.18; as = 1.03 ± 0.88 (B2) 

The per degree of freedom for the best-fit values of the parame- 
ters is found to be 1.05, which shows the fit is reasonably good. The 
difference between the data points and the fit is shown in Figure 
Bl. We can then represent the Itmiinosity distance obtained from 
the data by the function 

Qflt(2) = io»-2[-f.t(-)--^l (B3) 

Note that one needs to fix the value of M to obtain the function 
Qtit{z). For flat models, the value of M is well-constrained (the 
error being 1 per cent) - hence we can use the best-fit value of 
M = 23.95 to obtain 

_ 2(1 + 2.00^) 

- 1 + 1.032 ^^^^ 



